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Abstract 

Generalizing a recent convolution model of the nucleon which explicitly conserves 
charge and momentum we reexamine the influence of an enhanced pion field in nuclei 
on EMC and Drell-Yan ratios. Due to wave function renormalization constants the 
effect is more than 50% smaller than predicted by the standard pion excess model. 
In particular there is no discrepancy between the EMC data and our results in the 
x-region which was expected to be most sensitive to the pion. 



When in 1983 the EMC collaboration found an enhanced nucleon structure function 
F2{x) in iron at a; < 0.3 this was interpreted as a signature for an enhanced pion field in 
nuclei which was predicted already in the 70s . Starting from the pion distribution func- 
tion of the free nucleon, calculated in 1974 by Sullivan 0, Llewellyn Smith was able to 
explain the effect by the assumption that the iron nucleus contains about 0.15 additional 
pions per nucleon. Generalizing Sullivan's expression Ericson and Thomas |P showed that 
these extra pions result in a rather natural way from the attractive pion-nucleon interac- 
tion in a conventional meson-exchange description of nuclear matter. Consequently, after 
the enhancement of F2 could not be confirmed by later experiments [§], 0], it was con- 
cluded that there are no excess pions in nuclei, calling into question the meson-exchange 
picture itself @. These conclusions were corroborated by the non-existence of a nuclear 
dependence of Drell-Yan cross sections [|]. Recent attempts to solve the problem are 
based on a possible change of the gluon properties in nuclei p or the partial restoration 
of chiral symmetry with density [T^. In this communication we do not intend to argue in 



favor or against these suggestions. Instead we want to reinvestigate the question whether 
an enhanced pion field - if it exists - really leads to a detectable enhancement of the 
EMC and Drell-Yan ratios. In particular we want to discuss the role of wavefunction 
renormalization constants which have not been considered in the past. 

Starting point of our investigations is the two-phase model for the nucleon developed 
by Szczurek and Speth [|lT|. In this model the nucleon consists of a bare core and a 
meson-baryon phase. Here we want to restrict ourselves to pion and nucleon degrees of 
freedom which are supposed to be the main players in the region 0.1 < x < 0.3. Then 
the distribution function for a quark with flavor / in a proton reads: 

(1) 

For the sake of readability we do not write the Q^-dependence of the distribution functions 
explicitly. The two convolution integrals are due to the pion-nucleon phase or, more 
precisely, a vr+n phase with q = 2/3 and a vr^p phase with q = 1/3. The function f^^^ 



is the Sullivan expression for the momentum distribution of the pion while f^/^ is the 
analogous distribution function of the recoil nucleon. The sum is multiplied by an overall 
constant which normalizes the total probability of finding the dressed nucleon either 
in the bare nucleon phase or in the pion-nucleon phase to unity. It is given by 

Z-^^-l = Mf/^ = /' dy r/^iy) = f dy f'^{y) . (2) 

JO JO 

This is a consequence of the well-known number sum rules for valence quarks reflecting 
the conservation of baryon number, charge and strangeness: 

^ dx {q^j{x) - q^f{x)) = ^ dx {qlf,^^^{x) - qlf,^^^{x)) = N^f , (3) 

with A^^ = 2, A'^^ = 1 and Nf = 0. Note that eq. (2) also implies that the first moments 
of the distribution functions f^^^ and f^/^ have to be equal. 
Similarly momentum conservation 

Y,J^ dxxq^x) = dx X qlf,^^^{x) (4) 

leads to a relation between the second moments of /'^/^ and f^/^ [^]. In eq. (4) the 
flavor index / runs over quarks and antiquarks. In principle gluons have to be included, 
too. However, since experimentally the momentum fraction carried by quarks is about 
the same in a pion as it is in a nucleon (~ 47% at = 50{GeV/c)'^) we must demand 
momentum conservation for quarks and gluons separately when we want to keep the 
internal structure of the constituents unchanged. 

Higher moments of f^/^ are not related to f^^^ by conservation laws. However, 
according to the intuitive picture that the nucleon breaks up into a pion which carries 
the fraction y of the plus momentum of the original nucleon and a recoil nucleon which 
carries the missing fraction 1 — ?/ it is frequently assumed ||rT| 

fN/N^y^ = r/^(i-i/). (5) 

This automatically guarantees the validity of eqs. (2) - (4) but is not a necessary conse- 
quence of these relations. 

2 



Many features we want to discuss later on can be expressed in terms of the mean 
number of pions per nucleon, {71^^). For the free nucleon we get from eqs. (1) and (2): 

1 + Ml 

Note that even for large moments (ra^) is always less than one, in agreement with the 
two-phase picture of the nucleon we have started with: Since at any time there is at most 
one pion, (n^r) should never exceed unity. (Of course, for > 1, i.e. when the A^vr loop 
becomes more important than the bare nucleon phase the two-phase model itself becomes 
unrealistic and multiple pion phases have to be included. In practice, however, this does 
not happen.) 

The essential point of eq. (1) is that the renormalization constant Zn is an overall 
prefactor, i.e. not only the bare nucleon part but also the cloud contribution becomes 



renormalized. In contrast Melnitchouk and Thomas argue that one should multiply 
only the bare nucleon contribution by a normalization constant which is given by Zn = 
1 - Mf/^. This is of course what one gets when eq. (1) is expanded into a power series 
in the nNN coupling constant until order f?^^^. Numerically however, this makes a big 
difference. Note that in this case the mean number of pions per nucleon, (wtt) is identical 



with the first moment . The difference between ref. [0 and ref. |Tl| (i.e. our eq. (1)) 
is due to a different summation scheme of diagrams when the dressed nucleon propagator 
is calculated. The authors of ref. perform a systematic expansion in the coupling 
constant until quadratic order, thus taking into account only single one-loop diagrams. In 
contrast, the authors of ref. iterate the meson-baryon loops until infinite order. Very 



recently Holtmann showed by a comparison of the corresponding next order corrections 
that this is a significantly better approximation |jl3|. In any case, when we turn to the 
calculation of nuclear structure functions a strict expansion in the ttNN coupling constant 
is no longer a useful scheme. The medium modifications of the pion distribution function 
are mainly due to rescattering processes and Pauli blocking. So if we restrict ourselves 
to diagrams of order glj^j^ we should also drop all rescattering diagrams. In this case 
we would be left with Pauli blocking as the only medium effect which would lead to a 



depletion rather than an enhancement of the pion field. On the other hand, if we sum 
over an infinite set of rescattering diagrams as it is usually done we should also 
iterate the vrA^ loops which dress the nucleon. This means we should normalize like in 
eq. (1). 

In the pion excess model ||, [l^ it is assumed that the quark distribution functions of 



the bare nucleon and of the pion remain unchanged inside a nucleus. Thus in analogy to 
eq. (1) we write 

(7) 

with relations analogous to eqs. (2) - (4) in order to guarantee baryon number, charge and 
momentum conservation. Note that we still define the Bjorken variable x with respect to 
the nucleon mass. Therefore, unlike the free nucleon case, f^^"^ and /^M ^Iq not need to 
vanish for y >1 and in general a relation like eq. (5) can not be correct in a nucleus. 

On the other hand in microscopic models the upper limit ymax for non-vanishing 
does not really become equal to the nuclear mass number A as we have indicated in eq. (7) 



but remains rather close to unity. For the model presented in ref. for instance, we find 



Vmax = IH — 2- which is about 1.1 at nuclear matter density (The model of refs. [g, |16| even 



leads to Umax < 1 but that is an artefact of the nonrelativistic kinematics.). Furthermore 
the shape of the pion distribution function is not dramatically changed by nuclear effects 
and the main contribution still comes from the regime y < 1. Thus before we turn to a 
more realistic model we consider a schematic model where the pion distribution function 
in a nucleus is simply given by f"^^ multiplied by an enhancement factor A: 

/VA(^) = Ar/^(y). (8) 

For the distribution function of the recoil nucleons f^^^{y) we assume that a relation 
analogous to eq. (5) holds. 

Within this model we calculate EMC and Drell-Yan ratios for different values of A. 
For the irNN formfactor we choose a monopole parameterization with a cutoff parameter 



^ttNN = SOOMeV^. We assume the sea quark distributions of the bare nucleons and of 
the pions to be S'[/(3)-symmetric and the proton and the neutron as well as the different 
pions to be related by isospin rotations. We parameterize the quark distributions of the 
bare proton in the following simple form: 

^Qharei^) = Nj X^f (l - xff . (9) 

The parameters are chosen such that the experimental distribution functions of the dressed 
proton are roughly reproduced by eq. (1). We find af = 0.65 and f3f = 3 for valence 



up quarks, «/ = 0.65 and /?/ = 4 for valence down quarks and aj = 0, Pj = 9 and 
Nf = 0.17 for sea quarks. The normalization of the valence quarks is of course fixed by 
eq. (3) . For the pion we take the experimental quark distribution functions of ref . . 



The left panel of fig. 1 shows our results for the EMC ratio F^/F^ for isospin 
symmetric nuclear matter. For the deuteron nuclear effects are neglected. At small values 
of X the ratio becomes enhanced with increasing A whereas it becomes reduced at larger 
X. At X ^ 0.2 the EMC ratio remains unity independently of A. This behavior is easy to 
understand. Using F2{x) = we find for the EMC ratio at a; = 0: 

^^'(Q) = 1 + (in^) - (n^)) ^2^(0) no) 
Fi'iO) ^ ^^^^ ^ '^^^2^(0) ■ ^ ' 

This result is quite general and does not depend on the specific assumption eq. (8). For our 
schematic model we can also calculate the ratio at a; = 1. At this point the contributions 
of the pion and the recoil nucleon vanish and the EMC ratio is just the ratio of the 
renormalization constants: 

Fi'il) 1 - «) • ^ ^ 

Thus with an enhanced pion field ((n^) > {n^)) the EMC ratio is greater than one at 
a; = and less than one at a; = 1. So somewhere in between there must be a point Xq 
with F2{x,^ = F,^{xo). In our schematic model this point is given by by the condition 

Fi^K) = F^VK) (12) 



independently of A. We find Xo ~ 0.2. That means an enlianced pion field leads to an 
EMC ratio greater than one only for x < 0.2. In order to estimate whether this is a 
detectable effect we also plotted experimental data points in fig. 1. Note that at x < 0.1 
nuclear shadowing sets in which is not included in our model. So there is only a very 
small window 0.1 < x < 0.2 where an enhancement could be seen. For realistic parameters 
(A < 2) this seems to be very difficult. 

Here the wave function renormalization plays an important role. As we can see from 
eq. (10) the enhancement is directly related to the number of excess pions which is given 
by 

{n^) - {n^) = Zr,ZA{Ml"'-Ml"'). (13) 

Thus compared with calculations without normalization factors the effect is reduced by a 
factor ZnZa- Since Za < Zn for an enhanced pion field and Z^ — 0.7 for our parameters, 
the additional pion contribution is less than half of what it would be without normalization 
factors. It is also crucial that the renormalization includes the pion cloud, like in eq. (1). 
If we renormalized only the bare nucleon, like in ref. |]12|, we still would get eqs. (10) 
- (12). In this case, however, the pion numbers {n.,^) would be identical with the first 
moments and we would find eq. (13) without the factor ZaZ]^. 

In the right panel of fig. 1 we show the corresponding Drell-Yan ratios as a function 
of Xti the Bjorken variable for the target nucleus. The Bjorken variable for the projectile 
was kept fixed at Xp = 0.5. At this large value of Xp we are mainly sensitive to the target 
sea which is dominated by the pion contribution. Therefore this process is much more 
sensitive to detect an enhancement of the pion field although the effect is also strongly 
reduced by the normalization factors. 

After these schematic studies we want to make the model somewhat more realistic. 
The enhancement of the pion distribution function in nuclear matter has been calculated 
many times |^, |l^, [10|- In these models the bare pion propagator becomes dressed by 



nucleon hole and A hole excitations. As we have mentioned before, f^^^iy) does not 
have to vanish for y > 1 since a pion may carry more than the fraction 1/A of the plus 
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momentum of the nucleus. Therefore conservation of charge and momentum cannot be 
enforced by a relation like eq. (5). This can also be seen from more physical arguments: 
Since we are summing an infinite series of particle hole polarizations we should allow 
the virtual photon not only to couple to the recoil nucleon which corresponds to the 
first emitted pion but also to particle states which are created by absorbing a pion, i.e. 
to nuclcons which have interacted with another nucleon before. The most natural way 
to describe all these processes effectively is to incorporate them into Fermi motion and 
nuclear binding. 

We proceed as follows: Integrating eqs. (1) and (7) over x and using charge conser- 
vation, i.e. eq. (2) and the corresponding relation for the nuclear case, we can eliminate 
the recoil nucleon distributions as well as the distribution functions of the bare core. In 
this way we find expressions for / dxq^f'^{x) which only depend on the distribution func- 
tions of the (dressed) free nucleon and the pion distribution functions f^^^ and f^^"^. By 
construction the integrands already satisfy the requirements of charge conservation, but 
momentum is not yet conserved. In the next step we convolute the nucleonic part with 
a distribution function fpermii^) which describes the effects of Fermi motion and nuclear 
binding. This procedure does not change the first moments of the quark distribution 
functions but it changes the second moments as a function of a binding parameter 77. 
Thus by tuning rj to the appropriate value we can achieve momentum conservation while 
keeping the charge to be conserved. 

We find the following expressions for the quark distribution functions of a nuclear 
proton: 

Jx z z ^ Jx y y 

(14) 

with 

6 Jx Z Z Z 



7 



r^r/^(y)K"(-)-n-°(-))], 

y y y 

(15) 



A^(x) = —AP{x) and Aj(x) = in all other cases. 

Except for the normalization factors the two integrals in eq. (14) are identical with 
the standard pion excess model, as discussed e.g. in ref. [|l^. While the first integral 
describes the redistribution of the nucleons keeping the baryon number fixed, the second 
integral is the contribution from the excess pions. However, since the pion cloud around 
a proton carries a positive net charge, in order to conserve the total charge we have to 
correct for the extra charge of the excess pions: Due to the emission of additional vr"*' 
mesons the number of up quarks associated with the baryonic part of a nuclear proton is 
lowered while the number of down quarks is enhanced. This is the origin of Aj in eq. (14). 
Since this term does not contribute to the momentum balance, only its first moment is 
fixed and one can easily find other expressions which do the job as well as eq. (15). In 
practice, however, this does not really matter: Since for neutrons everything is just the 
other way around, in isospin symmetric nuclear matter the Aj-terms completely drop out. 
In nuclei with a small neutron excess they persist, but they are strongly suppressed. In 
fact, for ^^Fe we find that the influence of A/ can be almost neglected. Thus the main 



difference to earlier calculations g is the appearance of normalization factors in eq. 
(14) which reduce the pion enhancement in the same way as we have discussed for the 
schematic model. 

Our results for the EMC and Drell-Yan ratios {^^Fe compared with deuterium) are 
shown in fig. 2. The solid lines represent the ratios obtained from eq. (14) imposing 
momentum conservation. We take the experimental quark distributions of Owens for 
the pion as well as for the nucleon [|14|. The nuclear pion distribution function is 
calculated microscopically within the framework of a nuclear matter RPA calculation, 
including NN~'^ as well as AA^~^ excitations. Short-range correlations are taken into 
account by folding the potential with a correlation function gdr) = 1 — jo{<lc'^)- In the 
A^A^ channel we choose qc = 3.94/m~^ corresponding to the mass of the tu-meson. This 
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gives rise to a momentuiii dependent Migdal parameter with the hmit g'jsiN 

(g = 0) = 0.57. 

In the A^A and the AA channel qc is tuned to reproduce S'tvaIO) — q'aa^^) — V^, the 
classical Lorentz-Lorenz value. For details of the interaction see ref. [0, although in 
the present work we do not introduce density dependent masses and coupling constants 
which are discussed there. We use the nucleon distribution function fpermii^) derived 
from a Fermi gas model using the correct relativistic normalization [0. For the Fermi 
momentum we take kp = 260MeV which corresponds to p = 0.87po, the average density 
of the ^^Fe nucleus. The parameter t] is fixed by momentum conservation. In order to 
compensate for the additional momentum carried by the pions (~ 4.5%) we have to choose 
r] = 0.94. 

For comparison we also show the result obtained within the schematic model we dis- 
cussed above (dotted line). We choose A = 1.75 which corresponds to the enhancement of 
the pion field we find in the microscopic calculation. Of course, the schematic model is too 
simplistic to produce the rise of the EMC ratio seen at x > 0.7. However, for x < 0.6 and 
in particular for x < 0.3 - our region of interest - the two models agree quite well. This is 
also true more or less for the Drell-Yan ratio. Moreover, in a more detailed investigation 
(not shown in fig. 2) we find that a major part of the differences comes about from the 
fact that we do not exactly start from the same quark distribution functions of the free 
nucleon. Remember, in the schematic model we calculated q^{x) from eq. (1) using the 
simple parameterization eq. (9), whereas we directly used the experimental parameteri- 



zations [jT4| in the more realistic model. This is also the reason why the ratios do not 
exactly agree at x = where they are fixed by eq. (10). So with some more effort in 
finding a good parameterization for the quark distribution functions of the bare nucleon 
the dotted line would come even closer to the solid curve. 

The comparison shows that for this type of models the EMC ratio at small and in- 
termediate X is to a very large extent determined by the number of excess pions and the 
amount of momentum transferred from the nucleons to the pions: The first fixes the ratio 
at X = while the latter determines the slope of the drop with increasing x. Details 
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which are related to higher moments of the nucleon and the pion distribution functions 
are less important. Therefore we believe that the major conclusions we have drawn from 
the schematic model remain valid for more sophisticated descriptions. In particular it will 
be very difficult to produce any significant enhancement in the region 0.1 < 0.3. 

As mentioned above, in order to conserve the total momentum we have to choose the 
binding parameter to be 77 = 0.94. This corresponds to an average separation energy 
of (1 — rfjmN — 56MeV which is of course unrealistically large. In order to check the 
sensitivity to this parameter we give up the requirement of momentum conservation and 
choose rj = 0.97, the value proposed by Li et al. [l7|. The results are shown as the dashed 
dotted curves in fig. 2. Since we have not changed the number of excess pions the ratios at 
X = remain unchanged, but the EMC ratio drops off more slowly. However, for x < 0.4, 
i.e. in the region where the pion enhancement was supposed to show up, there is still no 
disagreement with the data. 

The momentum conserving EMC result (solid line) looks very similar to fig. 4 of ref. 
T6| . Note however, that our short-range repulsion is much weaker and more realistic |1^ . 



Without wave function renormalization this would lead to a much larger number of excess 
pions (as discussed in eq. (13)) and consequently to a larger ratio at small x. In order to 
demonstrate this we replace and Za in eq. (14) by unity. The results are indicated by 
the long and the short-dashed lines in fig. 2. When we enforce momentum conservation 
by choosing 77 = 0.88 we get the results indicated by the short-dashed line. Similar to 
the behavior discussed in fig. 1 the EMC ratio drops off very fast and becomes unity at 
X ~ 0.2. When we choose the realistic parameter rj = 0.97 we arrive at the long-dashed 
curve. This is the only case where the EMC-data are significantly overestimated in the 
whole area below x < 0.3. 

To conclude, we find that EMC experiments are not sufficiently sensitive to detect an 
effect due to a possibly enhanced pion field in nuclei. Within a schematic model we show 
that the enhancement of the pion distribution function must be unrealistically large to 
produce a significant deviation from the data in the region 0.1 < x < 0.3. These findings 
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are confirmed by more realistic calculations. Our results are based on a generalization of 
the two-phase model for the free nucleon of Szczurek and Speth In this model the 
pion contribution is strongly reduced by a renormalization constant. Other prescriptions 
which renormalize the bare nucleon part only are based on a power expansion in the 
coupling constant and seem to be unsuitable for nuclear matter. 

In contrast to EMC experiments, Drell-Yan experiments are sensitive to conclude that 
there is indeed no enhancement of the pion field in nuclei. Although the renormalization 
constants reduce the discrepancy by more than 50% the data are still overestimated in a 



standard RPA calculation. This leaves room for more unconventional explanations |10]. 

So far we have restricted ourselves to nucleon and pion degrees of freedom. Although 
heavier mesons are known to be more important in the intermediate regime x ~ 0.5, 
through the normalization factors and the momentum balance they also have some in- 
fiuence on the lower x region. We do not expect, however, that this will spoil our main 
conclusions. 
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and A. Szczurek for very useful information about their model of the nucleon. This work 
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thank J. Speth for his hospitality. 
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Figure Captions 



Fig. 1 EMC ratio (left panel) and Drell-Yan ratio (right panel) for isospin symmetric nu- 
clear matter calculated within the schematic model (eq. (8)) with A = 1 (solid line), 
A = 1.5 (dashed-dotted line), A = 2 (dashed line) and A = 4 (dotted line). The 
data points |^, are plotted in order to give an estimate for the experimental 
sensitivity. 

Fig. 2 EMC ratio (left panel) and Drell-Yan ratio (right panel) for ^^Fe calculated in a 
realistic model using eq. (14). The solid line shows the result for the binding pa- 
rameter 7] = 0.94, obtained from momentum conservation, while the dashed-dotted 
line is calculated with the more realistic value rj = 0.97. The dashed lines indi- 
cate the results we get in a calculation without renormalizing the pion contribution 
{Za = Zn = 1), short-dashed line: rj = 0.88, long-dashed line: r] = 0.97. For com- 
parison we also show the result obtained within the schematic model with A = 1.75 
(dotted line). The EMC data were taken from refs. (^°Ca) and (^^Fe), the 
Drell-Yan data were taken from ref . {^^Fe) . 
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